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Abstract For many quantum information implementations 
with trapped ions, effective shuttling operations are impor- 
tant. Here we discuss the efficient separation and recom- 
bination of ions in surface ion trap geometries. The maxi- 
mum speed of separation and recombination of trapped ions 
for adiabatic shuttling operations depends on the secular fre- 
quencies the trapped ion experiences in the process. Higher 
secular frequencies during the transportation processes can 
be achieved by optimising trap geometries. We show how 
two different arrangements of segmented static potential elec- 
trodes in surface ion traps can be optimised for fast ion sep- 
aration or recombination processes. We also solve the equa- 
tions of motion for the ion dynamics during the separation 
process and illustrate important considerations that need to 
be taken into account to make the process adiabatic. 



1 Introduction 

Significant progress has been made in quantum information 
processing with trapped ions [1,2 ,3 4], including entangle- 
ment gates [5 6 1, teleportation [7 . 8 .9 1 and quantum simula- 
tion ifToirrnfni. 

However it is not easily feasible to manipulate many 
qubits in a single trapping region. It would be useful if qubits 
can be stored in separate trapping regions (memory zones) 
and only be brought together in a single trap {processor zone) 
when quantum operations are required [ 13 14U15lfl6l . Shut- 
tling within an array of ion traps has been demonstrated suc- 
cessfully in linear arrays and through junctions [17,18,19 20 
|2TI . Two ions have also been reordered by rotation within a 
linear trap section [22]. Separation of two pairs of ions hav- 
ing different masses was also demonstrated ll23l . Ions are 
separated, re-combined and transported across the different 
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zones of an ion trap array by means of time-varying poten- 
tials on control electrodes. How to optimise electrode geome- 
tries for efficient ion separation and recombination has been 
discussed by Home and Steane [24] in general. House [25| 
studied surface-electrode ion traps analytically. Hucul et al. 
[26 1 and Reichle et al. Il27l have discussed the energy gain 
of trapped ions due to ion transport. We extend their findings 
to the ion separation process in surface trap arrays, providing 
a detailed description how such arrays need to be optimised 
to allow for efficient separation. Furthermore, we analyse the 
dynamics of the separation process by solving the equations 
of motion and present a description about the considerations 
that need to be taken in order to make the process adiabatic. 
There are two types of trap geometries, asymmetric ion traps 
B28U29IF301F3TI . where the electrodes lie in a plane and the ions 
are trapped above that plane and symmetric ion traps 11321 
[33U34lfr71l20l . where the electrodes are symmetrically posi- 
tioned around the position of the trapped ion. It is important 
to scale these architectures to trap and shuttle hundreds of 
ions for any useful quantum information processing to oc- 
cur. This article focusses on optimal geometries for asym- 
metric ion traps. Modern micro fabrication techniques are a 
promising approach to build such scalable ion trap arrays in 
which ions will be brought together and separated many times 
in processor zones to perform the gate operations. We will 
show that this is best attainable when the trap features are 
designed at the scale of the ion-electrode distance. In this ar- 
ticle we discuss the optimisation for the particular case of 
surface ion trap arrays. The speed of the adiabatic shuttling 
operation can be enhanced by maximising the secular fre- 
quency during separation and recombination inside the trap 
array [26]. The secular frequency depends on the applied 
voltages on the electrodes, the geometry of the traps and it 
typically increases for smaller ion-electrode distance. How- 
ever, at smaller scales, motional heating of ions becomes sig- 
nificant due to anomalous heating [35.36.37]. Cryogenic op- 
eration of ion trap chips may allow for small ion-electrode 
spacings as it is known to significantly suppress anomalous 
motional heating B36II37I . 
In surface ion traps the ion-electrode distance (ion height) de- 



Altaf H. Nizamani and Winfried K. Hensinger 




static potential provides confinement in the z-direction. The 
effective potential in all three directions is given by 13911261 . 



Figure 1 Illustration of the pseudopotential created above the sur- 
face of the trap electrodes when an rf voltage is applied on the rf 
electrodes while keeping the other electrodes at rf ground. x is 
the horizontal position of the ion with respect to the central ground 
electrode, h is the height of a trapped ion above the central ground 
electrode, and the escape point (turning point) for the ion is located 
beyond h. b and c are the widths of the rf electrodes and a is the 
separation between the rf electrodes. 



pends on the size and configuration of the electrodes M28II25I . 
The average ion life-time in a trap depends on the trap depth. 
One of the main challenges in surface traps is to achieve 
higher trap depths at larger ion-electrode distance, since such 
traps typically offer depths of about 1% that of multi layer 
symmetric traps of comparable dimensions [ 38 ] and the mag- 
nitude of the voltage that can be applied is limited by the ac- 
tual fabrication process. 

In this article we discuss how to design a surface trap array 
in which ions can be trapped at a maximum trap depth for a 
given ion-electrode distance and can be brought together and 
separated rapidly by adjusting static voltages on electrodes 
while maintaining the highest possible secular frequencies. 
The trap depth and secular frequencies are dependent on the 
applied voltages and geometric factors of surface trap geome- 
tries. As the applied voltages are limited by power dissipation 
and breakdown voltage of the trap electrodes, it is important 
to optimise the trap depth and secular frequencies by adjust- 
ing the dimensions of the electrodes. In Sec. [2] we show that 
the trap depth may be optimised at a given ion height by ad- 
justing the size and configuration of the electrodes. In Sec. 
[3j we show how to maximise the secular frequency during 
the separation and recombination shuttling processes by ad- 
justing the widths of the static potential electrodes for ion 
transportation in general, and fast ion separation processes in 
particular. In Sec. [4] we discuss the dynamics of the separa- 
tion process. We discuss constraints in the design of realistic 
trap arrays. We then compare two fundamental designs and 
present a guide to accomplish fast and adiabatic ion separa- 
tion. 



2 Optimisation of trap depth 

In a typical three-dimensional rf Paul trap the rf field provides 
trapping in the x and y dimensions (transverse axes) and a 
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where, e and m are the charge and mass of the ion be- 
ing trapped, r f (%) is the instantaneous electric rf potential 
when V r f = 1 V and Ot (x) is the static electric potential pro- 
duced by the ith static potential electrode when Vi = 1 V and 
X is the position vector. V r f is the peak rf voltage applied on 
the rf electrodes with drive frequency fi r f and the coefficient 
Vi(t) is the time varying voltage applied on the ith control 
electrode. The first part of Eq. Q] represents a pseudopoten- 
tial which can be created in a trap by applying an rf voltage 
on the rf electrodes while keeping the other electrodes at rf 
ground. In a surface trap, the position of the minimum of the 
pseudopotential or rf node where the ions can be trapped, is 
located at a distance h (ion height) in the y-direction and x in 
the ^-direction, as shown in Fig.Q] The escape point or turn- 
ing point of the pseudopotential shown in Fig. [T] is located 
above the trapping position. The position of the ion and the 
turning point can be found by calculating where the gradient 
of the pseudopotential is zero. In absence of any static poten- 
tial, the trap depth defined as the amount of energy needed 
for an ion to escape, can be represented by the difference be- 
tween the pseudopotential at the rf node (typically zero) and 
the turning point. In Fig.Q] the trap electrode dimensions are 
labelled as widths of the rf electrodes b and c and separation 
between the rf electrodes a. 

Small gaps between the trap electrodes in realistic geometries 
usually have negligible effects on trap parameters B40II251 . 
Therefore, the basis functions for the trap electrodes can be 
calculated using the analytical model incorporating the gap- 
less plane approximation [25 1. As House |25| suggested, if 
the origin of the coordinate system is located between the 
ground and left rf electrode as shown in Fig. Q] and the seg- 
mented static potential electrodes are considered infinitely 
long (in the x-direction), the rfnode can be found to be posi- 
tioned at, x = ac/(b + c) and h = y/abc(a + b + c)/(b + c) 
11251 . The ion should be at a reasonably large distance from 
the trap electrodes to reduce the effect of anomalous heating 
of the ion |37 36| and provide good laser access. The later 
requirement may be alleviated via the introduction of slots in 
the substrate that allow for optical access [42 .43 1. First we 
show how to maximise the trap depth for a given ion height 
h. 

Building on the discussion given by House 11251 . we re- 
express the trap depth S in terms of given ion height h and 
the geometric factor k for a given ion mass m and rf voltage 

Vrf, 
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where k is described as 
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Figure 2 Illustration of the principle axis rotation, (a) Orientation 
of the principle axis when equal width rf electrodes (b = c) are used, 
(b) The principal axis is rotated by using unequal width rf electrodes 
(b 7^ c). The horizontal shift Xd in the trapping position is caused 
by the configuration of the rf electrodes. To keep the static poten- 
tial electrodes symmetric around the trapping position in the axial 
direction, an extra ground electrode of width d is inserted between 
the narrower rf electrode and the static potential electrodes. 



Eq. (0 shows that the trap depth for a given ion height 
h, can be maximised by optimising the geometric factor k 
which is defined in Eq. ([3J. By choosing the appropriate elec- 
trode widths, the trap geometries can be optimised to achieve 
the maximum trap depth for a given ion height. Efficient laser 
cooling of an ion along all three directions of motion can only 
be achieved if the laser wave vector k has a vector component 
along all three principal axes. Therefore, in most ion trap ex- 
periments, rotation of the principal axes is required for effec- 
tive laser cooling of an ion in all three directions |44|. One 
of the techniques to achieve rotation of the principal axes is 
to use asymmetric rf electrodes which have different widths 
B31l as shown in Fig. [2] Unequal widths (6 ^ c) of the rf 
electrodes cause planar asymmetry in the ir-axis and set off 
the nonuniformity in the confinement field when equal static 
voltages are applied on opposite static potential electrodes. 
This issue can be managed by introducing an extra ground 
electrode of width d between the narrower rf electrode and 
the segmented electrodes parallel to the rf electrode as shown 
in Fig. [2] The width d of the ground electrode (in the gap- 
less approximation) may be chosen as the difference in the 
widths of the rf electrodes {Aw = \b — c\) plus the shift in the 
horizontal position Xd of the rf node, caused by the unequal 
widths of the rf electrodes. The width d can then be calculated 



as 



Ax 
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House 1 25 1 optimised trap depth for a given rf electrode 
separation a. However, we believe it is more useful to derive 
values for a given ion height since ion height is a major con- 
straint for many experiments due the occurrence of anoma- 
lous heating. The ratio C, = b/a between the rf electrodes 
widths, b and the separation between rf electrodes, a is use- 



ful to characterise k at a given ion height h. Using the ratio 
£, the geometric factor k can be parameterised for equal and 
unequal width rf electrodes as following 



C 2 (1 + 2C) 



4(i + C) 2 (i + C + %/T+2C) 2 



4C 2 (2 + 3C) 
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when c = b/2 
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Fig. [3] shows the trap depth geometric factor k as function of 
(. The solid curve in Fig. [3] shows k at a given ion height 
as a function of ( when b = c. The maximum of k can be 
found at £ « 3.68 for rf electrodes of equal width and for 
unequal rf electrodes when c = 6/2 at £ « 4.9 as shown in 
the dashed curve in Fig. [3] For the optimised values of £, the 
ion height above the electrodes is given by h w 1.43a for the 
equal width rf electrodes and h f=a 1.27a when c = b/2. The 
aim of an optimum trap design is to achieve the maximum 
trap depth at a given distance above the electrodes. In general, 
for optimised traps, the maximum trap depth decreases with 
increasing ion height h and scales approximately as ~ h~ 2 . It 
is important to note that anomalous heating of a trapped ion 
is proportional to ~ h~ 4 [35,36]. The trap depth increases 
with a decrease in ion-electrode distance, but the heating rate 
also increases with a decrease in the distance. Therefore, the 
aim of an optimum trap design is to achieve the maximum 
trap depth at a given distance above the electrodes. 



3 Optimisation of fast ion separation process 



In the ion separation process, initially, ions are trapped in a 
single potential well with secular frequencies ui x , ui y and cu z , 
normally (w x , u y ) > cj z , where uj x and uj v predominantly de- 
pend on the pseudopotential provided by the rf electrodes and 
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Figure 3 For a given ion height, the trap depth geometric factor k 
can be maximised for equal width rf electrodes (b — c) when the 
ratio between rf width and separation £ ~ 3.68 (solid curve). In the 
case when the rf electrodes widths are unequal and (c = 6/2), k is 
maximised when the ratio is £ « 4.9 (dashed curve). 



Altaf H. Nizamani and Winfried K. Hensinger 



12 



o 
0- 




-1.5 



Figure 4 Plots of the potential in the z-direction when a > 0, a — 
and a < 0. The potential wedge is created when a < (solid 
curve). 



lu z depends on the voltages applied to the static potential elec- 
trodes. The static voltages can be applied in such a way that a 
wedge potential can be created between the trapped ions and 
the single potential well can be pulled apart into two distinct 
potential wells or a double well in the z-direction of the trap 
geometry. The aim of an effective separation is that the ions 
remain trapped and acquire minimal kinetic energy during the 
separation process. Furthermore, we will show that, for adi- 
abatic separation to be possible, the separation process must 
be performed on the time scale of the minimum secular fre- 
quency during the separation. Following the theoretical work 
on ion separation by Home and Steane [24], the confinement 
potential near the centre of a trap in the z-direction can be 
analysed using a Taylor expansion as 



&{z,t) « 2ea(t)z 2 + 2e/3(t)z 4 



(6) 



Fig.|4]shows the plots of Eq. ((6]) when a > 0, a = and 
a < 0. A potential wedge can be seen in the middle of the 
trap when a becomes negative. Secular frequencies, u z , in 
the z-direction at each instance can be calculated using the 
equations derived by Home et al. Il24l 
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During the separation process, when the voltages Vi (t) on 
the static potential electrodes are varied in time, the quadru- 
ple term a crosses zero and at that point the secular frequency 
in z-direction, ui z , is at its minimum (w m i n ). At this point, the 
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Figure 5 Surface trap geometries with (a) outer segmented elec- 
trodes and (b) centre segmented electrodes. Unequal rf electrodes 
are used to provide rotation of the principal axes. 



The value of u>mm (when a — > 0) also sets an upper limit 
on the speed of the separation process. One should therefore 
aim to maximise uj m i n for a faster separation of ions. When 
a — >• 0, Wmin is only due to the contribution of quartic term (3 
in Eq. ©. 

Therefore summing up the conclusions from Home and 
Steane [24], in order for ion separation to occur, a trap de- 
sign must provide a negative value of a, when appropriate 
voltages are applied to the trap electrodes. Furthermore, the 
better trap design is the one which provides higher values for 
the quartic term /3 which in turn provides a higher value for 
UJmm during the separation process allowing for faster speed 
of the adiabatic shuttling process. From Eq. © one can see 
that, both a and /3 terms depend on applied voltages on the 
static potential electrodes and their dimensions. As we will 
discuss in Sec.|4] applied voltages on the electrodes are con- 
strained by the power dissipation and the breakdown voltage, 
therefore, the effective maximum value of f3 needs to be max- 
imised by optimising the trap geometry. 
We investigate two arrangements of static potential elec- 
trodes in surface traps for an effective and fast ion separa- 
tion process. The arrangement of rf and static potential elec- 
trodes in the two surface ion trap geometries are shown in 
Fig. |5] In Fig. [2a), the outer rf-ground electrodes are seg- 
mented to provide confinement in the z-direction. The ion - 
static potential electrode distance to the set of static potential 
electrodes on each side is made equal by inclusion of an ad- 
ditional axial ground electrode. Fig.[5jb) shows an alternative 
geometry where the central rf-ground electrode is segmented 
to provide axial confinement. The segmented electrodes are 
labelled as endcap, control and wedge. In both designs, ions 
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Figure 6 /? is plotted against the ratio of wedge electrode to endcap 
electrode width (W/E) (a) for the outer segmented electrode design 
and (b) for the central segmented electrode design. In both cases an 
optimal value of W/E is approximately 1.1. 
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Figure 7 /3 is plotted against electrode widths (W=C=E) in units 
of a (a) for the outer segmented electrode design, (b) for the central 
segmented electrode design. 



can be trapped and separated above the central electrode(s) 
by adjusting the voltages on the segmented electrodes [45 

[MET]. 

A surface ion trap lacks one of the reflection symmetries, the 
symmetry in the direction normal to the surface. Therefore, 
any applied voltage on the endcap electrodes and wedge elec- 
trodes during separation process can easily alter the position 
of a trapped ions above the surface [28 1 and push the ions 
out of the rf node position. The solution for this problem is 
to apply a negative voltage on two (in case of the design in 
Fig- 12b)) or more (in case of the design in Fig. [3 a)) control 
electrodes symmetrically around the rfnode. The voltage on 
the control electrodes can be maintained in such a way that 
the trapped ions always remain in the rfnode position, during 
the separation and shuttling processes. 

As we know from the previous discussion, an efficient 
separation process is dependent on the value of j3. Therefore 
the optimum trap geometry is one which provides maximum 
values for f3 and allows for negative values for a for certain 
applied voltages. As the values of the a and f3 parameters are 
also limited by the voltage constraints of the trap electrodes, 
in our numerical simulation we keep the applied voltage con- 
stant and equal to 1 V for endcap and wedge electrodes and 
— IV for control electrodes. This approach allows us to de- 



termine the dependency of the a and /3 terms on geometric 
factors of a surface trap only. The first step is to determine the 
optimum ratio (W/E) of wedge to control electrode. For this 
purpose we fix the width of endcap and control electrodes to 
equal widths and vary the width of wedge electrodes to max- 
imise j3. In Fig. [6] f3 is plotted against the ratio of wedge to 
endcap electrode widths for both trap designs. In these plots 
we can see that j3 reaches a maximum when the width of all 
the static potential electrodes are approximately the same. 
The widths of static potential electrodes should be chosen in 
such a way that they provide significant curvature of the po- 
tential at the centre of the trap in the z-direction and meet 
the conditions for effective ion separation by providing max- 
imum /3 and negative a. In Fig. [7] j3 is plotted as a function 
of electrode widths in units of rf electrode separation a for 
both surface trap designs. We can easily deduce from the plot 
shown in Fig. [3 a) that j3 can be maximised when width W 
of the segmented electrodes is sa 3.66a for the design shown 
in Fig. [2 a). Whilst, in Fig.[3b), j3 is maximum when W ~ a 
for the trap design shown in Fig. |5jb). The relationship be- 
tween (3 and the segmented electrode width in Fig. [7]for both 
designs also shows that a relatively larger width of electrodes 
is required to control the ion motion when the ion-static po- 
tential electrode distance is large. The ion - static potential 
electrode distance, for the outer segmented electrode geom- 
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Figure 8 (a) a and (b) /3 for outer and centre segmented geometries 
are plotted vs a. A unit voltage is applied on all trap electrodes. 
Comparison shows that the a and /3 values are higher in the centrally 
segmented trap geometry where ion electrode distance is small. 



etry, is ss \/[(6 + x )] 2 + [h] 2 , while for the central seg- 
mented static potential electrodes, this distance is « h. 
Furthermore, in Fig. [8] we compare both trap geometries fol- 
ia! and /?, where both terms are plotted as a function of rf 
electrode separation a (which dictates the height of the trap- 
ping position above the surface). The calculations are made 
with all other parameters optimised as explained above. By 
comparing both designs we find that the values for a and j3 
are approximately two order of magnitude higher for the de- 
sign shown in Fig.|3b). We can also observe a sharp rise in 
the values of a and /3 when the rf electrode separation de- 
creases below 50 fim, but the cost for this achievement can 
be severe because anomalous heating also increases rapidly 
when the distance between the electrodes and ion decreases. 



4 Ion dynamics problem during separation process 

Ion dynamics 

The aim of an optimal separation protocol is to allow the ions 
to be separated and transported to arbitrary locations within 
the trap array, whilst ensuring that the motional state of the 
ion does not change significantely after a shuttling operation. 
Hucul et al. |26| and Reichle et al. [27| have identified those 
constraints that ensure adiabatic transport of ions. Both have 



highlighted the importance of the inertial forcing of the ions 
at the beginning and end of a shuttling protocol. Hucul et al. 
showed the benefits of the hyperbolic tangent profile 



P h {t) =tanh[iV 



— ] 



(8) 



and Reichle et al. 11271 suggested an error function profile 

t-T 1 



P e (t)=Erf[n^=-] 



T 



(9) 



to shuttle the ions. In both equations (|8] and |5J, t is the in- 
stantaneous shuttling time and T is the total shuttling time. 
Functions with larger N or n produce a more gradual change 
for values of t close to and T, while incorporating larger 
rates of change for values of t close to T/2. If the N and 
n-parameters are selected appropriately, both profiles suf- 
ficiently resemble each other. Hyperbolic tangent functions 
take significantly less computation time than the error func- 
tion. Therefore we use the hyperbolic tangent time profile to 
analyse the dynamics of the ion separation processes. 
An arbitrary time dependent potential can be built using 
the basis functions for individual electrodes. The force on a 
charged particle can be calculated using the classical equa- 
tions of motion Il26l 



J2mx j + V j &{ Xj ,t) = 



(10) 



where the pseudopotential &(xj>t) is defined in Eq. Q] It is 
also possible to calculate the classical trajectories of ion mo- 
tion by solving Eq. ( TTOb numerically. High accuracy solutions 
of Eq. ( [Tol l provide the ion dynamics in the x, y and z direc- 
tions as a function of time t, which can be used to calculate 
the kinetic energy gained by the ion. In order to obtain the 
ion dynamics, we use a package called "NDSolve" provided 
in Mathematica-7 to solve these differential equations. 

Average motional energy 

In a quantum harmonic oscillator of frequency ui , the aver- 
age energy (E\ of level (n) is given by (E) = Tud((n) + h). 
In analogy to a classical harmonic oscillator and assuming the 
total energy of the ion is only due to its kinetic energy (which 
is maximum at the bottom of the potential well), the average 
motional quanta (n) for a trapped ion can be calculated as 



(«) s = 



2 mv t 
flWt 



(ID 



where m is the mass of the ion, v t is the maximum veloc- 
ity and uit is the instantaneous secular frequency. The kinetic 
energy of the ion in the frame of the pseudopotential well is 
due to its secular motion. By plotting the kinetic energy of the 
ion versus the shuttling time, the maximum kinetic energy of 
the ion at the start and end of the shuttling can be obtained. 
Hence, the change in the average motional state of the ion is 



<»>.= 



Final K.Err 
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hoj t 



(12) 
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where, K.E max is the maximum kinetic energy of the ion in a 
potential well. 



Motional heating caused by anomalous heating 

In a trap design with small ion-electrode distance, significant 
motional excitation of an ion can be caused by anomalous 
heating during the shuttling process. It was observed from 
experimental data that the heating rate (h) an is related to 
the ion-electrode distance scaling as h~ A and to the secular 
frequency scaling as cj~ 2 [36|. In order to provide a conser- 
vative approximation we can utilise a measurement we have 
recently carried out |46| and the scaling laws as stated above, 
to provide a realistic estimate of the motional state excitation 
of the Yb + due to anomalous heating 



(n) 



1.97 ±0.15 x 10 2t >m 4 Hz 



(13) 



where h is in units of micrometres and cj is in units of s _1 . 
We note that this expression is only valid for the particular ion 
trap and ion species it was measured for, however, it provides 
a reasonable estimate for illustration purposes. We note heat- 
ing rates can be substantially reduced by operation in a cryo- 
genic environment 1 36 , 37 1 as well as optimisation of elec- 
trode surfaces and materials. 

The motional quanta (n) gained from the anomalous 
heating during the shuttling process can be calculated by in- 
tegrating (nj over the shuttling time, 



(n) 

\ i an 



(n) dt 

\ I an 



(14) 



where, t and tf are the start and the end time for a shuttling 
process. This integral also takes into count the variation of 
the secular frequency ui z during the shuttling process. 
Therefore, the total number of motional quanta gained during 
ion transport is given by 

(n) = (n) s + (n) an (15) 



Separation in realistic trap geometries 

In order to demonstrate the importance of optimal trap ge- 
ometries for ion separation and recombination, it is useful to 
analyse the actual dynamics of the process. While the con- 
clusions to be obtained in this section are applicable for ge- 
ometries featuring a wide range of ion - electrode distances, 
we carry out actual simulations for a given set of example pa- 
rameters. Below we motivate the particular choice taken and 
note that the actual ideal set of parameters should be chosen 
depending on the particular fabrication process and other con- 
siderations such as whether the ion trap is operated in a cryo- 
genic environment, what motional heating rates are accept- 
able for the particular experiment, what ion species is used 
and what secular frequencies and trap depths are required. 



In Sec. |2] and [3] we investigated how ion trap electrode di- 
mension ratios can be optimised to maximise the k and (3 
parameters which provide for maximum trap depth and sec- 
ular frequency during separation. In order to carry out sim- 
ulations of the ion dynamics during separation we need to 
determine first a realistic set of voltages that can be applied 
to the chip which will determine the actual values of k and 
f3 in the shuttling process when using optimised geometries. 
The trap depth S, the rf trap stability factor q [39] and the 
power dissipation Pd in the trap all depend on the applied 
rf voltage V r f, the driving frequency fi r f and the given ion 
mass and the ion-electrode distance |48|. If there is no static 
voltage offset on the rf electrodes, the trap depth S and the rf 
trap stability q can be related as 



Yii 

2tt 2 



nq 



(16) 



where q — 2eV r / / (m,!7 2 f h 2 ) [48 1 and k is defined in Eq.[5] 
Microfabricated ion traps typically are limited in the amount 
of voltage that can be applied due to voltage breakdown via 
insulator bulk and surfaces. In order to achieve a deep trap 
one should therefore choose q as large as possible while still 
remaining safely inside the region of stability in parameter 
space. Utilising q sa 0.7 seems therefore a reasonable choice. 
Power dissipation within the ion trap can lead to heating of 
the trap chip, outgassing from trap material and eventual de- 
struction of the chip. Power dissipation Pd can be estimated 
as H47) 

P d « 0.5V? f O? f C 2 R (17) 

where, C and R are the trap capacitance and resistance 
respectively. For typical trap chip configurations a power 
dissipation of 3 W should not lead to a large temperature 
change of the ion chip. Considering a typical microfabri- 
cated chip with electrodes of thickness ~ 15 fim made of 
electroplated gold on a commercially available insulator 
wafer made of Silicon Oxide (SiC<2) (a similar technique 
is used by Seidelin et al. 11291 ). we can estimate typical 
values for the capacitance and resistance in such traps as 
R fa 0.5 Q and C w 20 pF, respectively. Setting q s=s 0.7 
for an 171 Yb + ion, power dissipation Pd < 3 W, assuming 
Q r f ~ 2-7T x 55 MHz, one could apply V r f of approximately 
450 — 500 V. Microfabricated ion traps typically feature a 
particular breakdown voltage (eg. largest voltage difference 
between adjacent rf and static potential electrodes). For this 
discussion we will assume a maximum voltage difference 
between adjacent electrodes on the order of 500 V and we 
will choose voltages applied to the electrodes accordingly. 
For the purpose of choosing an illustrative and realistic ion 
trap geometry we choose an ion-electrode distance suffi- 
ciently large to feature reasonably low motional excitation 
due to anomalous heating while allowing for reasonable trap 
depth and secular frequencies when realistic voltages are 
applied to the trap. If we choose an ion height of m 85 fim, 
a trap depth S ~ 0.32 eV and radial secular frequencies 
(uj x and uiy) of up to ~ 4.2 MHz for an 171 Yb + ion can be 
achieved by using the parameters mentioned above. Having 
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Figure 9 Voltage variation as a function of time for the control and 
wedge electrodes to produce the ion separation process in the outer 
segmented trap. 



chosen the height of the ion above the surface, the optimi- 
sation considerations in Sec. [2] and [3] uniquely determine all 
other electrode dimensions. The optimum rf electrode widths 
assuming unequal rf electrodes widths (b = e/2), chosen 
for principal axis rotation, are therefore b « 300 /im and 
c sa 150 ^m separated by a ground electrode of width 50 /im 
(i.e. a = 60 /im including the gaps of 5 /im). The optimum 
width of the control electrodes for the ion separation process 
are then W = C = E sw 220 /im for the design shown in 
|5j a) and W = C = E ps 60 /im for the design shown in 

Kb). 

Using these parameters we solve the equations of motion to 
calculate the ion dynamics and resulting overall motional ex- 
citation. While the actual results correspond to the particular 
trap parameters stated above, the conclusions obtained are 
applicable for all surface ion trap arrays. For simplification 
we will refer to the centrally segmented design (Fig. |5ja)) 
and the outer segmented design (Fig. |3b)) with constraints 
and dimensions explained above as Centre Segmented Trap 
and Outer Segmented Trap, respectively. 



Ion separation in the outer segmented trap 

First we discuss the dynamics in the outer segmented trap 
(Fig. 13 a)). We confine 1 Yb + ions using V r f ~ 450 result- 
ing in radial secular frequencies w x k, u) y ss 4 MHz. Axial 
confinement along the ^-direction is obtained by applying the 
static voltage of approximately +30 V on the endcap elec- 
trodes and approximately -34 V on the wedge electrodes re- 
sulting in a secular frequency in the z-direction of uj z ss 440 
kHz. Using these static voltages maximises the axial secu- 
lar frequency while still retaining sufficient trap depth. The 
ions are separated by changing the voltage on the wedge elec- 
trodes to approximately +50 V and the voltage on the control 
electrodes to approximately -48 V monotonically. At the end 
of the separation process, the ions are located in two distinct 
potential wells approximately 2 x W apart. The static volt- 
ages are chosen to achieve maximum secular frequencies dur- 
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Figure 10 Variation of the secular frequency ui z during the ion sep- 
aration process in the outer segmented trap using the hyperbolic tan- 
gent time profile with N = 3.0 (dashed) and N = 4.0 (solid). 
The secular frequency at the start and the end of the process is 
oj z /2n « 500 kHz. The lowest secular frequency experienced by 
the ions during the separation process is oj m i n /2n ~ 42 kHz. 



ing the separation process and to retain sufficient trap depth 
(at least 0.2 eV) during the transport. 

We use the hyperbolic tangent time profile for changing the 
voltage on the static potential electrodes in the ion separation 
process. To illustrate the functionality of the JV-parameter in 
the separation process, we use N — 3.0 and N = 4.0. Fig. 
[10] shows the voltages applied to the electrodes as a function 
of time. 

As discussed in Sec. [3] the secular frequency in the z- 
direction varies during the separation process. The variation 
of the secular frequency uj z for an 171 Yb + ion during the 
separation process followed by hyperbolic tangent time pro- 
file with N = 3.0 and N = 4.0 is plotted in Fig. [10] 
We can see that the secular frequency reaches a minimum 
w m ,„ w 2-7T x 42 kHz when the double well is about to appear. 
As shown in Fig.QJJ the secular frequency of the ions varies 
rapidly and at one point is at its lowest. Therefore, to reduce 
the energy gain during the separation process, the ions should 
be separated slower than the time scale of the minimum sec- 
ular frequency 27r/u; m i n . The gain in motional quanta (n) is 
plotted in Fig. [TTJas a function of shuttling time for hyper- 
bolic tangent with N = 4.0 and N = 3.0. The solid curves 
represent the motional quanta gain due to the shuttling pro- 
cess (n) , whilst the dashed curves show the motional quanta 
gain caused by anomalous heating (n) , Note we use a par- 
ticular heating scaling law (eq. [T31 which is only valid for a 
particular ion trap and ytterbium ions. We only use this law 
for illustration purposes, the use of other ions species and trap 
materials will result in different absolute values of motional 
excitation even though the observed trends will remain the 
same. The crossing points between these two curves provide 
a reasonable idea of the minimum motional excitation that 
can be achieved. While longer shuttling times will reduce the 
motional excitation that results from the actual shuttling pro- 
cess, they will increase motional excitation due to anomalous 
heating. Therefore finding the crossing point between the two 
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curves provides for the optimal shuttling time scale. We stress 
that the actual minimum achievable excitation in a particular 
shuttling process is dependent on the actual motional heating 
rate, the figures here only serve to illustrate the principle. 
The gain of the motional quanta (n) also depends on the 
TV-parameter of the hyperbolic tangent profile. This can be 
explained by the profile of the variation of the axial secular 
frequency during the separation process with TV = 3.0 (solid 
curve) and TV = 4.0 (dashed curve) shown in Fig. [10] The 
graph shows that the ions spend a relatively long time in the 
lower frequency region during the separation process when 
the value of the TV-parameters is smaller therefore being sub- 
ject to more motional excitation via anomalous heating. Fig. 
[TTI shows that approximately the same number of quanta is 
added at the crossing points for the separation profiles with 
TV = 3.0 and TV = 4.0, but the duration of the separation 
is smaller in case of TV = 3.0. Optimising for the best value 
of TV allows for small gains in achievable lowest motional 
excitation and shuttling speed. 



Ion separation in the centre segmented trap 

Next we discuss the centre segmented trap illustrated in Fig. 
12b). Ytterbium ions are initially stored in a single potential 
well applying V r f « 500 V and a static voltage of approxi- 
mately 8 V on the endcap electrodes. The static voltages are 
chosen to confine the ions with approximately maximum sec- 
ular frequency while not making the overall potential anti- 
trapping and retaining at least 0.2 eV overall trap depth. With 
these voltages applied, the maximum secular frequency in ax- 
ial direction is lo z s=s 1 MHz and the radial secular frequen- 
cies are lj x s=a uj y w 4.3 MHz. We separate the ions by ad- 
justing the voltage on the wedge electrode to approximately 4 
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Figure 11 Gain in the average motional quanta (ro) by the ion in the 
External Segmented Trap after the separation process as a function 
of shuttling time. Hyperbolic tangent time profiles using TV = 4.0 
and TV = 3.0 are used to change the voltage on the control elec- 
trodes. The solid lines represent the best fit to the average motional 
quanta (71) resulting from the shuttling process and the dashed 
lines shows the gain of (n) from motional heating in the trap. 
The crossing points set lower limits for the total gain in (n) during 
the shuttling process. 
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Figure 12 Voltage variation as a function of time for the control 
and wedge electrodes to produce the ion separation process for the 
Centre Segmented Trap. 



V and the voltage on the control electrodes to approximately 
-3.8 V. In order to optimise the TV-parameter, we carry out 
simulations using the hyperbolic tangent time profile with 
TV = (3.0, 4.0, 4.5). Similarly as in the case for the exter- 
nal segmented geometry, the secular frequency in z-direction 
varies during the separation process as shown in Fig. Qj] The 
lowest secular frequency U) m i n for an Yb + ion during the sep- 
aration process is approximately 230 kHz. Fig. [T2l shows the 
voltages applied to the electrodes as a function of time. 
The average motional quanta (n) gained by the ion after the 
separation process are plotted as a function of total shut- 
tling time in Fig. [14] The solid curves represent the mo- 
tional quanta gain due to the shuttling process (n) , whilst 
the dashed curves show the motional quanta gain caused by 
anomalous heating (rij . The crossing points of (n) (solid 
lines) and the (n) (dashed line) set the lower limit of the 
total average motional quanta (n) gained by the ion during 
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Figure 13 Variation of the secular frequency cj z , during the ion sep- 
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Figure 14 Gain in the average motional quanta (n) by the ion in the 
Centre Segmented Trap after the separation process as a function of 
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the gain of (n) from motional heating in the trap. The crossing 
points set lower limits for the total gain in (n) during the shuttling 
process. 



the separation process. There are only minor differences for 
the different N parameters. 



Comparison 

By analysing the actual dynamics of the ion separation pro- 
cess using realistic examples we can learn a lot about optimal 
separation. While one may assume optimal ion trap geome- 
tries are only useful to provide for faster adiabatic separation 
processes, our results show that optimal geometries may in 
fact be a prerequisite for adiabatic separation due to the exis- 
tence of anomalous heating caused by fluctuating charges on 
the trap electrodes. Furthermore, we show that the minimum 
gain of total motional quanta (n) during the separation pro- 
cess is much lower in the centre segmented electrode geom- 
etry due to the higher achievable values of the axial secular 
frequency uj z during the separation process. In fact, the centre 
segmented geometry allows for much faster separation with 
overall motional excitation still remaining very small. It is 
also possible to achieve higher secular frequencies in the cen- 
tre segmented geometry while still retaining sufficient overall 
trap depth. 

These results also demonstrate the importance of optimisa- 
tion of electrode dimensions as derived in Sec.|2]and[3] A ge- 
ometry with optimised trap depth provides the capability for 
applying larger static voltages (before the overall potential 
becomes anti-trapping), due to the deeper trap depth which 
in turn provides for larger axial secular frequencies and faster 
ion separation. Particularly in the case of the outer segmented 
trap geometry, it is important to use optimised electrode di- 
mensions in order to partially compensate for the in-principle 
limitations caused by anomalous heating in order to estab- 
lish at least near adiabatic operation. We stress that the actual 



values of estimated total motional excitation only serve illus- 
trative purposes and are expected to significantly vary when 
using different ion species and ion traps. 



5 Conclusion 

We have demonstrated that effective and fast ion separation 
in scalable surface ion traps at maximum trap depth can be 
achieved by optimising sizes and arrangement of the ion trap 
electrodes within a surface ion trap array. We have calculated 
the optimum ratio of the widths of rf electrodes over their 
separation for the maximum trap depth at a given ion height. 
The trap parameters a and (3 which characterise the secular 
frequencies during the separation process can be maximised 
by optimisation of the electrode dimensions. We have solved 
the equations of motion for the dynamics of ion separation 
and illustrated the importance of optimised electrode config- 
urations. A separation process performed with higher secu- 
lar frequencies adds less amount of energy to the ions. We 
have shown that centrally segmented ion trap geometries are 
superior in their performance compared to outer segmented 
geometries. Centrally segmented geometries allow for signif- 
icantly smaller overall motional excitation and also provide 
for much faster adiabatic separation processes. In fact, de- 
pending on the actual experimental conditions, they may even 
be a prerequisite to accomplish adiabatic separation. Due to 
the much simpler fabrication of outer segmented geometries, 
these may nevertheless be a geometry of choice. In that case 
our article illustrates both the importance and design such a 
geometry with optimal trap dimensions. 
Ion trap arrays are of significant importance for the imple- 
mentation of scalable quantum technology with trapped ions. 
Ion separation within such arrays may likely play a critical 
role and our paper shows how this process can be accom- 
plished optimally. 
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